Viscoelastic Taylor-Couette instability of shear banded flow 
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We study numerically shear banded flow in planar and curved Couette geometries. Our aim 
is to explain two recent observations in shear banding systems of roll cells stacked in the vorticity 
direction, associated with an undulation of the interface between the bands. Depending on the degree 
of cell curvature and on the material's constitutive properties, we find either (i) an instability of 
the interface between the bands driven by a jump in second normal stress across it; or (ii) a bulk 
viscoelastic Taylor Couette instability in the high shear band driven by a large first normal stress 
within it. Both lead to roll cells and interfacial undulations, but with a different signature in each 
case. Our work thereby suggests a different origin for the roll cells in each of the recent experiments. 



The Taylor-Couette instability that arises when a New- 
tonian fluid is sheared between concentric cylinders has 
a long history in classical hydrodynamics, dating back 
to the groundbreaking paper of G. I. Taylor in 1923 [l[. 
The effect is inertial in origin: fluid is forced centrifugally 
outward along the radial direction r, and recirculates via 
roll cells stacked in the vorticity direction z. In contrast, 
inertia is usually negligible in the non-Newtonian flows of 
complex fluids. However the last two decades have seen 
considerable interest in an inertialess, viscoelastic Taylor- 
Couette (VTC) instability 0, originating instead in the 
hoop stresses (first normal stresses) that arise when a 
polymeric fluid is sheared in a curved geometry. These 
squeeze fluid radially inwards, again triggering an insta- 
bility that leads to roll cells stacked along z. 

Another intensely studied flow phenomenon in complex 
fluids is that of "shear banding" [3j , in which an initially 
homogeneous shear flow gives way to a state of coexist- 
ing bands of unequal viscosities and internal structuring, 
with layer normals in the flow-gradient direction r. Close 
analogies exist between this non-equilibrium transition 
and conventional equilibrium phase coexistence, both ki- 
netically and in the ultimate banded state. There are also 
fundamental differences, e.g., in the way the coexistence 
state is selected in the absence of a free energy minimisa- 
tion principle Indeed, beyond the basic observation 
of banding, an accumulating body of data shows that 
many (perhaps most) shear banded flows show compli- 
cated spatio-temporal patterns and dynamics [Bj]. These 
are often associated with roll cells stacked along z, in 
both curved Couette and planar Q flow geometries, 
the origin of which remains unclear. 

In this Letter, we give the first theoretical evidence 
to suggest that, in a curved flow, these rolls can arise 
via a mechanism in which the high shear band, once 
formed, develops large enough first normal stresses to 
trigger a further bulk patterning instability of the VTC 
type within itself. We further show that this bulk in- 
stability disappears below a critical value of the cell cur- 
vature q, depending in a quantifiable way on the fluid's 
constitutive properties, and consistent with the known 
criterion Q that VTC instability occurs above a criti- 



cal value of 75 1 / 2 . At small curvatures, however, a dif- 
ferent instability emerges, driven by a jump in second 
normal stress across the interface between the bands 0] . 
This also leads to roll cells, but with identifiably different 
properties from those of the bulk instability. By mapping 
in this way a phase diagram containing these two sepa- 
rate instabilities - a bulk VTC-like instability of the high 
shear band, and an instability of the interface between 
the bands - we suggest different origins for the different 
experimental observations of roll cells in Refs. 0, 0|- 

Our study thereby brings together three different hy- 
drodynamic instabilities in complex fluids: shear banding 
itself; instability of an interface between bands; and, for 
the first time theoretically, a bulk VTC-like instability of 
one band. We hope thereby to stimulate further experi- 
ments, in a family of flow cells of different curvatures, to 
verify (or otherwise) our findings with regards this strik- 
ingly rich array of hydrodynamic instabilities. 

Governing equations — The generalised Navier-Stokes 
equation for a viscoelastic material in a Newtonian sol- 
vent of viscosity r\ and density p is 

p(d t + v.V)v = V.(T - PI) = V.(S + 2r/D - PI), (1) 

where v(r,t) is the velocity field and D is the symmetric 
part of the velocity gradient tensor, (\/v) a p = d a vp. We 
assume zero Reynolds' number p = 0; and incompress- 
ibility d a v a = 0. The quantity £(r,i) is the viscoelastic 
contribution to the total stress T(r,t), assumed to obey 
diffusive Johnson-Segalman (DJS) dynamics 0] 



(<9 t + t> • V) £ = a(D-S + S-£>) + (S-n + n-S) 

(2) 



-2GD + — V 2 £. 

r r 



Here a is a slip parameter, which must lie in the range 
— 1 < a < 1; G is a plateau modulus; r is the viscoelas- 
tic relaxation time; and $1 is the antisymmetric part of 
the velocity gradient tensor. The diffusive term V 2 £ is 
needed to correctly capture the structure of the inter- 
face between the bands, with a slightly diffusive interfa- 
cial thickness 0(1); and to ensure unique selection of the 
shear stress at which banding occurs [4[ . 
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FIG. 1: Bulk flow curve of ID base state for q — 0.0 (left), q = 
0.16 (right). Solid symbols (i), (ii), (I), (II) are for reference 
later in the manuscript. Styles of vertical lines indicating low 
and high shear bands are used consistently with Fig. [2] 



Geometry and boundary conditions — We study flow 
between concentric cylinders of radii r = {i?i,i?2}, in 
cylindrical coordinates r, 9, z. The inner cylinder rotates 
at speed V; the outer is fixed. We denote the cell cur- 
vature by q = (Ri — R^j/Rx, the flat limit of planar 
Couette flow corresponding to i?i — > oo at fixed R2 — R1, 
and so to q — > 0. The natural bulk rheological variables 
are then the cylinder velocity V, which is the main ex- 
perimental control parameter; and the torque T = q 2 T, 
which we have rescaled by q 2 to remain finite and equal 
to the shear stress T xy as q — ► 0. We assume invari- 
ance with respect to 9 and study the model's dynamics 
in the flow-gradient/vorticity plane r—z, which is the one 
most commonly imaged experimentally. At the cylinders 
we assume boundary conditions of zero flux normal to 
the wall n.VE Q) 3 = V a, (3 for the viscoelastic stress, 
although others are possible [10]; and no slip or perme- 
ation for the velocity. The height of the simulated region 
in the vorticity direction is L z , with periodic boundaries. 
We use units in which G = 1, r = 1 and R2 — R\ = 1. In 
these we choose parameter values 77 = 0.05, as suggested 
experimentally @, 0] ; a cell height L z = 2.0; and an order 
of magnitude estimate I — O(10~ 3 ) [111 ]. 

Numerics — We change variables from (r, z) to (p, z) 
where p = ]n(r/Ri)/ ln(i?2/-Ri), so mapping our curved 
geometry onto an effectively flat one, amenable to sim- 
ulation on a regular rectangular grid in which the only 
effect of curvature q is to generate extra source terms in 
the governing equations. Our study holds for arbitrary 
values of q. Components of the governing equations are 
then extracted, discretized on a grid (pj, Zi), and evolved 
in time as described for the special case of planar Couette 
flow, q = 0, in Ref. 0. Results below have timestep and 
grids (Dt,Dp,Dz) = (0.0002, a/512, 1/512) with a = 1 
(resp. 1/2) for unbanded (resp. banded) base states, 
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FIG. 2: Shear rate and normal stresses in low and high shear 
bands vs. slip parameter a (flat geometry, q — 0.0). 



convergence checked against smaller (Dt, Dp, Dz). Be- 
fore generating new data we checked that our code repro- 
duces known results for (i) the instability of an interface 
between shear bands in planar Couette flow Q; (ii) ID 
banded states in curved Couette flow (iii) dispersion 
relations for the onset of VTC instability in the Oldroyd 
B model for small q 2]; (iv) neutral stability curves for 
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this Oldroyd B VTC instability for finite q 

Results: ID banded states — First we discuss the re- 
sults of one-dimensional (ID) calculations that artificially 
assume translational invariance in z, allowing structure 
only in the main shear banding direction r. 

In the flat limit q — ► 0, force balance dictates that the 
shear stress T xy is uniform across the gap. Within the 
assumption of a similarly homogeneous shear rate, the 
homogeneous constitutive relation is given by T xy (j) — 
^xy(i) + Vl where T, xy (^) follows from solving Eqn. [2] 
subject to invariance in time and space. See the thin 
solid line in Fig. Q] (left). For an applied shear rate in the 
region of negative slope, homogeneous flow is linearly un- 
stable and the system separates into coexisting bands at 
a selected stress T se i V 1 — a 2 — 0.483. The steady state 
bulk flow curve therefore shows a plateau in the banding 
regime (thick solid line in Fig.Q] left), at this single value 
of the stress for which a stationary interface can exist be- 
tween bands. For non-zero cell curvature the shear stress 
varies as T xy ~ 1/r 2 . This forces the high shear band to 
reside next to the inner cylinder. As this band expands 
outwards with increasing applied shear rate, the overall 
torque increases to ensure that the interface between the 
bands remains at the selected stress T se \. Consequently 
the stress "plateau" in the banding regime acquires a 
non-zero, curvature dependent slope: Fig.Q] right. 

In principle, these ID base states must be calculated 
separately for each set of values (i],a). However we fix 
1] = 0.05 throughout, and further exploit the fact that 
base states for all values of a collapse onto a single mas- 
ter scaling state when expressed in terms of the variables 
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FIG. 3: Left: Maximum growth rate for an applied shear 
rate just into the high shear branch vs. the first normal stress 
scaling variable 1/(1 — a), for two different cell curvatures: 
(qr,7(l - a 2 ) 1/2 ) = (0.16,8.58) and (0.0,6.87). Right: Same 
for an applied shear rate 7(1 — a 2 ) 1 / 2 = 1.91 in the banding 
regime. States above the horizontal line are unstable, uj* > 0. 
Value of 1/(1 — a) in Fig. [4] is shown by an arrow. 
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FIG. 4: Left: Maximum growth rate for a scaled applied 
shear rate 7(1 — a 2 ) 1//2 = 1.91 in the banded regime, vs. cell 
curvature q. The slow interfacial instability at q = 0.0 (inset 
shows zoom near origin) is eliminated by curvature. At higher 
curvature, the high shear band displays a VTC-like instability. 
Right: corresponding wavevector. Symbols (i), (I) match 
those in Fig. [JJ here for one fixed value of the first normal 
stress scaling variable 1/(1 — a) — 416 (arrow in Fig. [3}. 



7VI - a 2 , VVl - a 2 , rVl - a 2 , ^xyVl - a 2 , H xx (l — a), 
and + a). Accordingly, the scaled flow curves 

in Fig. Q] represent all values of a in the allowed range 
— 1 < a < 1. Reading off the scaled shear rates of each 
band (vertical lines in Fig. [TJ left) we can then easily 
extract their true shear rates as a function of a (Fig. [2J 
top) . The corresponding first and second normal stresses 
are also shown in Fig. [5] (for q — 0, but the trends are the 
same for a curved cell, q ^ 0). For a — ► 1, the shear rate 
and first normal stress become large in the high shear 
band. We therefore anticipate, and will show below, that 
for sufficiently large cell curvatures q this leads to a bulk 
instability of the VTC kind in this band. Note that we 
are using a as a parameter that de facto controls the 
strength of the normal stresses, which is in reality set 
by a combination of material properties such as worm- 
like micellar length and degree of branching. We thereby 
hope to stimulate new experiments on shear banding flu- 
ids with different constitutive properties. 

Results: interfacial vs. bulk instability — We now 
move beyond the assumption of ID, performing 2D sim- 
ulations in the (r, z) plane. In each run we take as our 
initial condition, or "basic state" , a ID solution as just 
discussed, to which we add a 2D perturbation of tiny 
amplitude. By monitoring the early time evolution of 
the Fourier modes exp(ikz) exp(wt) of this perturbation, 
we can extract the dispersion relation ui{k) that char- 
acterises the linear stability of the basic state. We then 
read off the maximum uj* of this function. A value uj* > 
(resp. uj* < 0) shows the basic state to be linearly un- 
stable (resp. stable) to 2D perturbations. 

We start by considering a basic state comprising ho- 
mogeneous flow on the high shear branch in a flat ge- 
ometry q — 0, denoted (ii) in Fig. [TJ The corresponding 
dispersion relation reveals this to be linearly stable, as 



expected: its peak uj* < (solid line in Fig. O left). In 
contrast in a curved device, denoted (II) in Fig.[TJ the ba- 
sic state becomes linearly unstable when the first normal 
stress exceeds a critical value, consistent with the onset 
of a bulk VTC-like instability (dashed line in Fig.H left). 
This presence of a VTC-like instability in flow states that 
reside fully on the high shear branch leads us to expect a 
bulk instability of the VTC kind in the high shear band 
of a banded flow, in a curved cell q > 0. Accordingly, 
we now turn to the stability properties of a shear banded 
basic state. We start with the flat case q = before 
turning to the main situation of interest with q > 0. 

Shear banded flow in a flat geometry q — is already 
known [8] to show an instability driven by a jump in 
second normal stress across the interface between the 
bands, leading to undulations along the interface with 
wavelength Awl. This is seen in Fig. [3] (solid line, right) 
for the base state (i) in Fig. [TJ The maximal growth rate 
uj* is weakly positive in the range < a < 1 [161] . 

So far, then, we have brought together for the first time 
in the same model three instabilities already documented 
separately in the literature: shear banding itself 0]; a 
VTC-like instability of a strongly sheared polymeric ma- 
terial with a large first normal stress in a curved flow 
cell Q; and the instability in a flat geometry of an in- 
terface between shear bands with respect to undulations 
with wavevector in the vorticity direction driven by 
the jump in second normal stress across the interface. 

Our most significant result, however, is to report the 
stability properties of a shear banded state in two situ- 
ations not previously studied theoretically - in a curved 
flow cell, and when the high shear band has large first 
normal stress. We do so to demonstrate a bulk VTC-like 
instability of the high shear band. 

As shown by the dotted and dashed lines in Fig. [3J 




FIG. 5: Greyscale snapshots of T, xx on the ultimate attractor 
in the p — z (horizontal- vertical) plane for a shear rate in the 
banding regime jyl — a 2 = 1.91. Left: interfacial instability 
in a flat cell q = with 1/(1 — a) = 1.43. Others left to right: 
bulk instability of the VTC kind in the high shear band for 
a large value of the first normal stress scaling variable 1/(1 — 
a) = 416, as in Fig. H for curvatures q = 0.115, 0.13, 0.16. 



right, curvature suppresses the interfacial instability just 
discussed, dragging the weakly positive mode for q = 
below the axis uj* = in the left hand part of the plot. 
For large first normal stresses, though, we find a cross 
over (interchange of dispersion maxima) to a different 
instability: now of the bulk VTC kind in the high shear 
band. For a fixed value of the first normal stress scal- 
ing variable 1/(1 — a) the trend with cell curvature is 
shown in Fig. [4] left. Again the interfacial instability (i) 
is suppressed by non-zero curvature q > 0, before the 
bulk VTC instability onsets in the high shear band at 
q i=s 0.112. At the interchange, the wavelength of the 
fastest growing mode switches (Fig. [4] right) from Awl, 
consistent with an interfacial instability, to A ~ 1/8, con- 
sistent Q with a VTC instability in a bulk flow phase of 
width i=s 0.3, for this applied shear rate. Greyscale snap- 
shots on the ultimate nonlinear attractor for each kind 
of instability are shown in Fig. [5] The stress signals of 
the rightmost of these appear chaotic, consistent with the 
observation of complex roll cell dynamics in Ref . [6| . We 
defer to future work a detailed study of the temporal roll 
cell dynamics in our simulations. 

Summary and outlook — For large first normal stresses 
and cell curvature we have demonstrated a bulk instabil- 
ity of the VTC kind in a high shear band. The phase 
boundary (Fig. [6]) scales as q ~ (1 — a) ~ -Vf 1 ~ 7~ 2 , 
consistent with the criterion for VTC instability known 
for bulk unbanded flow [2|. In contrast, for small curva- 
tures and first normal stresses we find an undulatory in- 
stability of the interface between the bands as reported in 
Ref. Q for q = 0. An important additional contribution 
of this work has been to show this interfacial instability 
to be suppressed by cell curvature, our data suggesting 
the same threshold scaling q ~ (1 — a) as for VTC. 

The aim of this Letter has been to shed light on two 
different recent experimental observations of roll cells in 
shear banded flows @, 0] ■ Nghe et al. Q demonstrated 
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FIG. 6: Phase diagram in the plane of curvature q vs. first 
normal stress scaling variable 1/(1 — a). Dashed lines show 
power —1 as a guide to the eye. 

an instability leading to rolls stacked along the vorticity 
direction z in a rectilinear microchannel. The lack of 
cell curvature suggests these experiments to correspond 
to (the pressure driven equivalent of fl3| ) our (boundary 
driven) calculations for q = in Fig. [51 and so to a linear 
interfacial instability. Indeed, it has long been known 
that flows with parallel streamlines are linearly stable 
with respect to bulk perturbations, although a nonlinear 
(subcritical) instability cannot be ruled out [fH ]. 

Lerouge and coworkers j|| demonstrated roll cells 
stacked along z in a Couette cell of curvature q — O(0. 1). 
They further estimated their high shear band to satisfy 
the criterion for bulk VTC instability, which would corre- 
spond to the top-right of our Fig.[6j Indeed, the complex 
dynamics of Ref. Q do perhaps suggest VTC, although 
the wavelength ps 1 seems more consistent with inter- 
facial instability. To resolve this issue fully, it would 
be extremely interesting to perform a series of experi- 
ments scanning right across Fig. [6l Most immediately, 
this could be achieved by using a family of flow cells 
of different q. While the stability gap in Fig. [5] appar- 
ently covers an impractical range of curvatures, we would 
expect smaller values of the high shear viscosity 77 (not 
uncommon experimentally but difficult to access numeri- 
cally) to narrow this gap significantly, such that it indeed 
has the potential to be spanned by a family of cells. 

Other open questions include the interaction these roll 
cells with unstable modes of wavevector in the flow di- 
rection [l5| ; and the effect of stick-slip dynamics at the 
wall on this rich array of hydrodynamic phenomena. 
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